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Abstract 

The horizontal component amplitudes observed by ground-based observatories of the INTER- 
MAGNET network have been used to analyze the global pattern variance of the solar diurnal 
variations. Data from magnetic stations present gaps in records and consequently we explored 
them via a time-frequency gapped wavelet algorithm. After computing the gapped wavelet trans- 
form, we performed wavelet cross-correlation analysis which was useful to isolate the period of 
the spectral components of the geomagnetic field in each of the selected magnetic stations and 
to correlate them as function of scale (period) with the low latitude Vassouras Observatory, Rio 
de Janeiro, Brazil, which is under the South Atlantic Magnetic Anomaly (SAMA) influence and 
should be used as a reference for an under-construction Brazilian network of magnetic obser- 
vatories. The results show that the records in magnetic stations have a latitudinal dependence 
affected by the season of year and by the level of solar activity. We have found a disparity on the 
latitudinal response at southern and northern hemispheres during solstice, which is expected due 
to the asymmetry of the Sq field. On the other hand, records from magnetic stations located at ap- 
proximately the same latitude but at different longitudes presented peculiar dissimilarities. These 
results suggest that quiet day patterns and the physical processes involved in their formation still 
remain an important issue, far from being fully understood. 

Keywords: Magnetogram data, H component, Quiet days, Gapped Wavelet analysis, Wavelet 
Cross-correlation. 



1. Introduction 

The daily variations o f the g eomagnetic field were discovered by the english researchers 



Graham and Watchmaker ( 1724) through the observation of a compass needle motions in 1722. 



Since then, it is well known that a typical spectrum of those magnetic variations is composed by a 
few harmonics of the 24-h period (12, 8 and 6-h). They are described as "quie t daily geomagnetic 
field variations" - referred to as "Sq" for "solar quiet" field Campbell ( 1989b . 



The traditional method of calculating the baseline for the quiet day variations is to use the 
five quietest day for each month for each magnetic observatory. In this work, we propose a new 
way to deal with quiet periods by eliminating the disturbed days using a multiscale process. To 
accomplish this task, we study the harmonics of the solar diurnal variations using hourly data of 
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the H component and explore it via a gapped wavelet technique (described in lFrick et all 1 19971) . 
The gapped wavelet technique is suitable for analysis of data with gaps. In our case, datasets, 
from various mag netic stations that present gaps in the records in different periods of time (see 



Frick et al.l (11998b for mathematical motivation and proofs) were analyzed. 



Originally wav elet transform s we re applied in geophysic s to an alyze seismic signals in the 
pioneer works of Morlet ( 19831) and iGrossmann and Morletl (119841) . Nowadays, the use of the 
wavelet technique has exponentially grown in many different areas. A similar trend is noted in 
the application of the wavelet cross-correlation techniqu e, which has been used by many work- 
ers. It was first in troduced bv lNesme-Ribes et al.l (11995b to study two time series of solar data. 



Frick et al.l (1200 II) . in space research, analyzed optical and radio maps of the nearby spiral galaxy 



NGC 6946 using the correlation between images as function of scale. 

Nevertheless, the wavele t cross-correlation can be used for many other different approaches 
and purposes. For instance. lOczeretko et al.l (12006b used the wavelet cross-correlation method to 
assess the synchronization between contractions in different topographic re gions of the uterus, 
as a p otential diagnostic alternative in several pathologies. More recently. iRehman and Sid diqi 
(2009) researched correlations between pairs of time series of meteorological parameters such as 
pressure, temperature, rainfall, relative humidity and wind speed to study the climatic dynamics 
at different locations. 

This work aims mainly to highlight and interpret the solar diurnal variations at a Brazilian 
station compared to other twelve magnetic stations reasonably well distributed over the whole 
Earth's surface. By applying gapped wavelet transforms to these signals, we were able to ana- 
lyze both the frequency content of each signal and the time dependence of that content. After 
computing the wavelet transform, we performed wavelet cross-correlation analysis, which was 
useful to isolate the period of geomagnetic spectral components in each station and to correlate 
them as function of scale (period). 

The rest of the paper is organized as follows: Section [2] is devoted to explain the principal 
mechanisms and aspects of the solar magnetic variations. In Section|3] the analyzed period and 
data are presented. Section [4] describes, divided in subsections, the used methodology: Sec- 
tion |4T|presents a brief description of continuous wavelet transforms, Section l4~2l is devoted to 
introduce the gapped wavelet analysis, and Section l4~3l to establish the wavelet cross-correlations 
and to explain how they can be quantified. In Section [5] the results are discussed and, finally, 
Section |6]brings the conclusions of this work. 



2. The Physics of the Solar Magnetic Variation 

The major driving force for quiet day field changes seem s to arise from the dynamo-current 
process in the ionospheric E region between 90 and 130 km dStewarn. il 882). Tidal winds move 
the ions across the Earth's magnetic field producing electro-magnetic forces (emfs). Those emfs 
drive electric currents in the conducting E region which give rise to daily variations in the mag- 
netic field measured at the ground level (for details see IChapman and Bartelsl, 119401) . Through 
these mechanisms, two vortices of currents are induced, one in the northern hemisphere (clock- 
wise) and another in the southern hemisphere (counterclockwise). At the same time, a strong 
eastward electric jet is formed throughout the equatorial region. 

There are three factors that affect the dynamo process: the ionospheric wind, the ionospheric 
conductivity and the geomagnetic field configuration. The wind and the co n ductivity vary sea 



sonally due to their dependence on the solar zenith angle (ICampbelll 119891) . IZhao et al.l d2008l) 
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concluded that the correlation between the Sq amplitude and solar zenith was better in high 
latitude than in low latitude regions due to the effect of the prenoon-postnoon asymmetry of Sq. 

The Sq field variation has a main spatial dependenc e on latitude and is affected b y other factors 
including epoch of the year and level of solar activity. IChapman and S tagg ( 1929) observed two 
kinds of regular changes in the solar diurnal magnetic variations on quiet days: annual variation, 
that affects both the type and the amplitude during each year ; and so l ar acti vity variation, that 
affects fundamentaly the amplitude along each sunspot cycle. iTakeda (1999) estimated that the 
intensity of the Sq currents in high solar activity was about twice larg er than in low so lar activity. 

Many other researchers studied the variations of the Sq, in cluding Iffibberd ( 19851) that exam- 
ined the annual, semi-annual and even the whole solar cycle. SteningT i 197ll) examined seasonal 
variations and longitudinal inequalities of the electrostatic- field in the ion osphere by looking at 
its electric conductivity and the Earth's main magnetic field. ITake da (2002) showed solar activity 
dependence of the Sq amplitude, and explained this effect through the ionospheric conductivity. 
Takedal (2002) also compared the amplitude of the Sq for the same value of conductivity. The 



seasonal variation is seemingly due to differences in neutral winds or to the magnetic effect of the 
field-aligned current (FAC) flowing between the two hemispheres generated by the asymmetry 
in the dynamo action. The FACs are controlled by interplanetar y magnetic fi elds (IMF) and its 
electric fields can directly penetrate to the equatorial ionosphere dSastril ll988). 

Some ev idence of the influence of oceanic tides on the magnetic daily variation has been 
obtained bv lLarsen and Coxl(l966). They found small semidiurnal variations of the Z component 
at a coastal site (Cambria, California) and at two island stations (Honolulu and San Miguel) that 
could not be explained by the atmospheric tidal theory. They suggested that these variations 
must be due predominantly to oceanic tides. It is important to mention here that the conductivity 
of the ocean does not vary significantly with time, unlike the ionospheric conductivity. As a 
consequence, the seasonal variation of the oc eanic contribution is expected to be smaller than the 



ionospheric con tribution dCueto et aU 12003). 



In his review, Price d 19691) suggested that the interpretation of field variations at a particular 
observatory must be done carefully. Although the ionosphere current system may have a fairly 
simple world-wide pattern, the relatively small variation of the conductivity distribution of the 
Earths surface will introduce a corresponding small variation into the distribution of the induced 
currents. 



3. Magnetic Data 

In this section, we first describe the data used to study the quiet day variations. We also discuss 
the considerations used in data treatment for the different magnetic stations. 

3.1. Dataset 

In this paper, we used ground magnetic measurements to study the correlation between Vas- 
souras (VSS - Brazil) and other 12 different magnetic stations, which we selected to provide a 
reasonably homogeneous distribution of measurements on the Earth's surface. The distribution 
of the magnetic stations, with their IAGA code, is given in Fig.Q] The corresponding codes and 
locations are given in Table [T] This work relied on data collections provided by the INTER- 
MAGNET programme (http://www.intermagnet.org). 

We chose the period of geomagnetically quiet days (through the Kp index, to evaluate mag- 
netospheric disturbances) to study seasonal and solar activity variation. The data interval used 
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Figure 1 : Geographical localization of the stations used in this work and their respective IAGA code. 
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in this work was from 1999 to 2007, almost the whole solar cycle 23 using the available dataset. 
In this data interval, the higher solar activity occurred between 2000 and 2001; and the lower 
solar activity, between 2005 and 2007. For seasonal variations, one year was divided into three 
seasons: June solstice, equinoxes (March and September) and December solstice. 

To determine the disturbed days from the quiet days the Kp index was used. The global Kp 
index is a number from to 9 obtained as the mean value of the disturbance levels within 3-h 



interval observed at 13 subauroral magnetic stations (see Bartels, 1957). Due to our interest in 
studying quiet days variations, we considered quiet days, only those days in which the sum of the 
Kp indices was less than or equal to 24 and the Kp index was not higher than 3. In our approach, 
the disturbed periods were eliminated and considered as gaps. Therefore, in principle only the 
magnetic effects of the lunar and solar contribution remained in the magnetic records. 



3.2. Data Treatment Considerations 



We used low- and mid-latitude magnetic stations so that the magnetograms for these days 
should be not seriously affected by auroral magnetic disturbance s and the Sq would correspond 
to a regularly recurring phenomena in each day in these stations (lPricel ll969). 

In this work, we focus on the solar diurnal variations, and particularly, its diurnal and semid- 
iurnal period spectral components. These features were greatly explored over the last decades, 
as well as their spatial depend ence on lati t ude and others factors, including epoch of the 
year and level of so l ar activity dPricel 1 19691: ISteningl 119711 Iffibberdl 119851 ICampbellL [l997 : 
Le Sager and Huangil2002l:lTakedall2002l) . 



is We took into consideration the two well known facts about the Sq field cited bv lPricd (1969): 



1. It is largely a local time field that can be roughly represented by a current fixed to the Sun, 
so we used the data at Universal Time (UT) to observe the global variance of the Sq field. 



2. It is influenced by the Earth's main field throughout the latitudes between 60"N and 60"S, 
so we did not use magnetic stations of high latitudes (auroral zone). 
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Table 1 : INTERMAGNET network of geomagnetic stations used in this study. 



Station 


Geografic coord. 


Geomaj 


'iietic coord. 


Lat.(°) 


Long.(") 


Lat.(°) 


Long.C) 


AMS 


-37.83 


77.56 


-46.07 


144.94 


ASP 


-23.76 


133.88 


-32.50 


-151.45 


BEL 


51.83 


20.80 


50.05 


105.18 


BMT 


40.30 


116.20 


30.22 


-172.55 


BOU 


40.13 


-105.23 


48.05 


-38.67 


CLF 


48.02 


2.26 


49.56 


85.72 


CMO 


64.87 


-147.86 


65.36 


-97.23 


EYR 


-43.42 


172.35 


-46.79 


-106.06 


HER 


-34.41 


19.23 


-33.89 


84.68 


HON 


21.32 


-158.00 


21.59 


-89.70 


KAK 


36.23 


140.18 


27.46 


-150.78 


SJG 


18.12 


-66.15 


27.93 


6.53 


vss 


-22.40 


-43.65 


-13.43 


27.06 



Source: http://wdc.kugi.kyoto-u.ac.jp/igrf/gggm/index.html (2010) 



At high and polar latitudes, the geomagnetic field lines are aligned almost vertically so the 
ionospheric currents are joined with the field-aligned currents (FAC), and the electro dynamics 
is do minated by the influence of the solar-wind-magnetosphere interaction processes dYomotol 
2006). On the other hand, the ionospheric current at mid dle and low latitude is generated by the 
influence of tidal winds, i.e., ionospheric wind dynamo ( Sastrii 1988 ). 

Our primary interest was to correlate the response of the geomagnetic field at Vassouras to the 
other twelve previously chosen magnetic stations. The Vassouras Magnetic Observatory, is lo- 
cated in Rio de Janeiro, Brazil, under the South Atlantic Magnetic Anomaly (SAMA) influence. 
It has been active since 1915 and is a member of the INTERMAGNET programme. One of the 
peculiarities of VSS is its location, at a low latitude, where the H component is essentially the 
same as the total geomagnetic field. In forthcoming years, a Brazilian network of magnetometers 
will be implemented and VSS could be used as reference. 

To fulfill our purpose, we used hourly mean value series of the H geomagnetic component. 
Some magnetic stations have available solely the X component, so we converted the X com- 
ponent of XYZ system to the H comp onent of the HDZ system of vector representation of the 
Earths magnetic field (as described in Campbell, 19891) . We used X — H cos(D), where X is 
the vertical component in the XYZ system, H is the horizontal magnitude and D is the angular 
direction of the horizontal component from the geographic north (declination). In principle, this 
system conversion does not affect our results because we only used magnetic stations of low- and 
mid-latitudes. 

The magnetic stations that use the XYZ system are AMS, ASP, BEL and CLF. The conversion 
of systems was performed in BOU, CMO, EYR, HON and SJG for the data between 1999 and 
2003 and in BMT and KAK for the data between 1999 and 2004. 

An adequate knowledge of the daily variation field and a full understanding of the associated 
phenomena can only come from exte nsive and de tailed analysis of the mean hourly values of the 
magnetic elements at many stations d Price , Il969h . So the availability of these data in the World 
Data Centers that collect and distribute these data is of enormous help. 
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4. Methodology 



151 The method used in this study is based on continuous wavelet transforms. In this section, we 

152 first introduce the concepts of the continuous wavelet transform (CWT) and its properties. Fol- 

153 lowing, we introduce the gapped wavelet analysis and emphasize its improvements over the 

154 CWT. After that, we describe how we extracted the information of the observed data using 

155 wavelet cross-correlation analysis. 



157 
158 



174 



4.1. Continuous Wavelet Transform 

The wavelet transform was introduced at the beginning of the 80s of t he XX century bv lMorlet 



( 1983 ), who used it to study seismic data. iGrossmann and Morlet ( 1984 ) improved the windowed 



159 Fourier transform and they constructed the continuous wavelet transform (CWT). The idea was 

160 to change the width of the window function accordingly to the frequency of the signal being 

161 considered. The CWT is an integral transform, and it yields an affine invariant time-frequency 

162 representation. A wavelet transform can measure the time evolution of the frequency transients 

163 within the signal. In this transform, we use as basis, a set of functions called wavelets. 

164 Formally, in order to be called wavelet, a function if/ musts satisfy some conditions. First, the 

165 area under the curve described by the function iff, must be zero, 



iff(f)dt = 0. (1) 



This condition is known as the admissibility condition and assures that the wavelet have a 
mean zero. Second, the function should be localized in both, physical and Fourier spaces (time 
and frequency), i.e., the Heisenberg relation At Alj = const must be satisfied. 

The CWT of a time series / is defined by the integral transform, 

W{a,b)= J f(t)r a ,„(t)dt, (2) 

where * represent the complex conjugate. This function W(a, b) represents the wavelet coeffi- 
cients that are function of both time and frequency (position b and scale a). It is obtained by 
dilating and translating the chosen analyzing wavelet and convolving it with the time serie (/(?))■ 
The function ifr is expressed as 

xj/ ab {t) = ^—\li\ t -—\, a>0, aeR, (3) 



where the pre-factor \a\ is introduced in order to guarantee that all the scaled functions \a\ 

175 have the same energy. The parameter b means a translation in time. However, the time-frequency 

176 resolution of the wavelet transform depends on scale a. We have good time localization (however, 

177 poor frequency localization) when a is small and vice-versa. 

178 There are many possible analyzing wavelet. The choice of these wavelet depends on both: data 

179 and analysis objectives. Therefore, does not exist the best wavelet for signal analysis in general. 
iao For instance the Haar wavelet has a compact support in physical space and a larg e support in 

181 Fouri er space, expected by the Heisenberg's uncertainty principle (as discussed in Daubechies, 

182 Il990h . For relatively long data lenghts, smooth wavelet without a compact support (Maar and 

183 Morlet) can be used for computing CWT in terms of the Nyquist sampling theorem. The Maar 

184 wavelet is well localized in both, time and frequency, and it has a better localization in time than 
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the standart Morlet wavelet. Morlet wavelet is a complex wavelet and produces complex wavelet 
coefficients that gives us amplitude and phase information. In summary, the characteristic of the 
wavelet transform depends on the chosen analyzing wavelet. For instance, to study a time series 
with abrupt variations (steps) the Haar wavelet may be the most convenient; but for smoother 



time s eries Morlet and Maar wavelet are more recommended (see for instance Domingu es et al 
2005, and references therein). 
We used the Morlet wavelet, 



i/r(f) = exp 



192 
193 
194 
195 
196 
197 



202 
203 
204 
205 

206 
207 
208 
209 
210 
211 
212 
213 
214 



YcT 1 



exp(/ wo t), 



(4) 



with i — y—l, (±>o is the frequency parameter and cr is the time resolution parameter. The 
adjustable parameter cr gives the optimal time-frequency resolution. Small values of cr give 
better time resolution, while large values improve frequency resolution. In the gapped wavelet 
case, the problem of the admissibility condition breaks down when cr is b elow 1 . It must b e 



avoided, and therefore, a wide range of cr should be used (as described in ISoon et all 119991) . 
Here, after testing, we adopted a>o — 6 and cr = 1 that also provides the better possible time- 
frequency equilibrium. In this case also, the Morlet wavelet is suited very well for experimental 
data analysis because it has a Gaussian envelope. It allows reach a reasonable compromise 
between the time-frequency resolutions. 

It is possible to analyze a signal in a time-scale plane. In the wavelet analysis it is called 
the wavelet scalogram. In analogy with the Fourier analysis, the square modulus of the wavelet 
coefficient is used to provide the energy distribution in the time-scale plane. As in the Fourier 
analysis, we can also explore the dominant periods (or frequecies) of the time series through the 
global wavelet spectrum. It helps to understand the behaviour of the energy at a certain scale. 

4.2. Gapped Wavelet Analysis 

Magnetograms are finite length observational data series and may contain gaps of various 
sizes. To reduce gap problems, we use the ga pped wavele t analy sis a technique introduced by 



£ga 



Frick et al.l (119971) and afterwards improved in lFrick et alj (119981) . The mathematical properties 
of this algorithm were studied in detail. In this transform, the admissibility condition is broken 
when the wavelet overlaps data gaps. The leading idea of the gapped technique is to restore the 
admissibility condition by repa iring in some way the wavelet itself. 

Following [Frick et al.l dl997l) . we separate the analyzing wavelet in two parts, the oscillatory 
part h(t) and the envelope <p(f), 



if,(t) = h(t) ip(t), 
h(t) = exp(/ too t), 



(5) 
(6) 



scm = c\p| — 



(7) 



When the wavelet is disturbed by the gap, we can restore the admissibility condition by in- 
cluding a function K(a, b) in the oscillatory part of the wavelet, 



if/(t , b, a) 



t-b 



K(a, b) 



t-b 



(8) 
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and requiring, 



222 
223 



\l/(t) dt = 



(9) 



The introduced function K(a, b) can be determined for each scale a and position b from (0 
and®. 

It was shown that this technique not only suppresses the noise caused by the gaps and bound- 
aries, but improv es the accuracy of frequency determination of short or strongly gapped signals 
dFrick etal.[|l998l) . 



4.3. Wavelet Cross-correlation Analysis 

The approach of this work is to use the wavelet cross-correlation to study the correla- 
tion between a pair of mag netic data from d ifferent stations as a function of scale (see 
Nesme-Ribes et al.l 11995) and Fric k et al. (2001) for more mathematical details): 



C(a) 



J <W\(a,t)<W\(a,i)dt 
(/ «Wi(fl, f) 2 dt J <W< 1 {a, t) 2 dtf 



(10) 



where 'Tt / ,(a, t) - \Wi(a, t)\ - \Wi(a, t)\, Wj are the wavelet coefficients and W, is the arithmetic 
mean in time for i = 1 or 2. We can also obtain the determination coefficient as function of the 
scale a. It is defined as the square of the correlation coefficient. One of the reasons to use the 
coefficient of determination instead of the correlation is to compute the statistics to determine 
the size or magnitude of the relation between two variables. It is interpreted as the percentage of 
variance in one variable that is predicted or explained by the other. The correlation coefficient 
permits the prediction of a score on one standardized variable from another, but this prediction 
will be meaningful only if the relation between the two variables is approximately linear (see for 



237 instance lOzerl 1 1 985l and references therein). Though in space geophysics both the determination 

238 and correlation coefficient are used. In our case, we preferred to use the determination coefficient 

239 due to this interpretation of variance. 

240 As we mentioned above, in the wavelet representation the scale (a) defines the distribution of 

241 energy of the whole signal through time. The wavelet cross-correlation allows us to check the 

242 interaction between two sets of data for each considerated scale. The scales are chosen in such 

243 a way that they make possible to characterize the dominating periods in the geomagnetic data 

244 spectrum. 

245 5. Results and Discussion 



246 In order to exemplify the wavelet scalogram expected from t he analysis of a daily magnetic 

247 variation, we defined a synthetic Sq similarly to the described in lMaslova et al.l d2010l) . 



248 and the noise R(t) as 



2nt 



Sq(t) = sin — — + sin — — 



24 



4nt 



R(t) = 0.5* V(f)*sin| — 



24 



2nt 



(11) 



(12) 



249 



The good Sq estimative Q(t) is given by 



Q(t) = [S q (f) + R{t)] *U(f), (13) 

250 where the day to day variability was introduced by the noise variables V(t), U (f) and R(t), where 

251 V(t) was uniformly distributed on [0, 1] and U(t) on [0.5, 1.5]. 

252 Fig. |2fa), from top to bottom, shows the synthetic geomagnetic signal, the wavelet square 

253 modulus (scalogram) and phase with Morlet wavelet. It was possible to identify the two most 

254 prevailing periods (12 and 24-h) due to the quiet variations. Also, we can identify the change of 

255 frequencies introduced by the noise, which indicates that the wavelet analysis is useful for this 

256 multi-scale phenomena. The phase change was also analyzed, which allowed us to retrieve the 

257 oscillatory behavior of the data. 

258 In Fig.|2b), we used the same synthetic signal used in Fig. [2 a) with addicional random gaps 

259 in the signal. The areas with gaps were signalized by arrows and dashed lines. We applied the 

260 gapped wavelet in this signal. In this new analysis, it was also possible to identify the quiet 

261 geomagnetic variations with periods of 12 and 24-h. 

262 In an attempt to fill the gaps, we used cubic splines which led to an underestimation of high 

263 frequencies. As the signals may have many spectral components, only by using gapped wavelet, 

264 the high-frequency part of the spectrum can be reconstructed. For this reason, the gapped wavelet 

265 has some advantages over interpolation. However, the gapped wavelet has some limitations, it 

266 can only be used if the gap is not larger than the period of analysis. 

267 After we performed the CWT and gapped wavelet in the synthetic signal, we concluded that 

268 the gapped wavelet fulfilled our purpose. Then, we started calculating the wavelet coefficients 

269 for each chosen station. 

270 Fig. [3]shows an example of CWT applied to a real signal, in this case, a magnetogram data. In 

271 this analysis, we used the H component obtained in June, 2007 at VSS station. Fig. |3ja) shows 

272 the H component, (b) the wavelet square modulus (scalogram), (c) the global wavelet spectrum 

273 (total energy in each scale) and (d) the phase. Areas of stronger wavelet power are shown in 

274 dark red on a plot of time (horizontally) and time scale (vertically). In this wavelet power, it is 

275 not necessary to include the cone of influence because the gapped wavelet technique not only 

276 suppresses the boundary an d gaps in bo t h high and low frequencies, but also estimates better the 

277 frequency of the signal (see lFrick et al.l d 1998) for mathematical proofs). Using gapped wavelet, 

278 we were able to determine both the dominant modes of variability and how these modes vary in 

279 time. In Fig.[3jb), the 24-h period shown in red is clearly visible and the 12-h period in the reddish 

280 color appears less dominantly compared to the 24-h. In Fig. [3c), the global wavelet spectrum 

281 shows the pronounced increase of the total energy in the 24-h period and a less pronounced but 

282 still noticeable in the 12-h period. 

283 The gapped wavelet used on the sythetic signal was then applied in real data as shown in 

284 Fig. |4] One might notice that Fig. [4] is very similar to Fig. [3] the difference is that now we 

285 introduced gaps due to disturbed days. As our interest was only to study the quiet day variations 

286 we considered the disturbed days (where the Kp > 3) as gaps. Even with the disturbed days, the 

287 global wavelet spectrum of Fig. |4] still shows a pronounced increase of energy in the 24-h and 

288 less dominant in the 12-h period. The inclusion of gaps in the data due to the removal of the 

289 disturbed days does not affect the analysis of the diurnal and semidiurnal variations. With the 

290 gapped wavelet coefficients in hand, we excluded the coefficients that corresponded to the gaps 

291 areas in the scalograms. The excluded areas are presented as white rectangles in Fig. [4] 

292 After the gapped wavelet analysis, we have calculated the cross-correlation functions C (a), for 
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293 each pair formed by Vassouras and one of the twelve chosen magnetic stations. This procedure 

294 was done for the years 1999 to 2007 during the equinoxes and solstices. These calculations were 

295 performed in order to understand the global response of the Sq variation at different locations. 

296 Fig. |5]shows the modulus of the correlations functions (for June, 2007) for the twelve magnetic 

297 stations. The correlation graphics are displayed from top to bottom and divided in two blocks, 

298 in the same order as in Table Q] These selected stations are divided in low latitude (ASP, BMT, 

299 HER, HON, KAK and SJG), medium latitude (AMS, BEL, BOU, CLF and EYR) and high 

300 latitude (CMO). In the vertical axis, we present the correlation coefficients and in the horizontal 

301 axis, the scale (period in hours). The periods of 24, 12, 8 and 6-h are highlighted with dashed 

302 lines in order to facilitate visual inspection and comparison (see Fig.|5J- It was possible to verify 

303 that for most of the pairs C (a) varies considerably with scale. We observed that the correlation 

304 between two geomagnetic data sets is scale dependent. They were strongly correlated or less 

305 correlated depending on the scale. We also observed that the correlation was usually larger for 

306 the first harmonics of the diurnal variations, and consequently got smaller and smaller for the 

307 following harmonics. Through Fig. [5] we determined the correlation for 24 and 12-h periods and 
3oa also calculated the determination coefficients. 

309 As mentioned before, the principal driven forces of the Sq currents are due to the daily thermal- 

310 tidal motions and the thermospheric winds moving the ionization through the Earth's magnetic 

311 field. During the solstices, there is a wind asymmetry between the Northern and Southern hemi- 

312 spheres due to the contrast in summer- winter radiation. In Fig. [5] we expected a small correlation 

313 between VSS and the magnetic stations located in the Northern hemisphere (BEL, BMT, BOU, 

314 CMO, HON, KAK and SJG) as a result of the month of June be a solstice: the thermosphere 

315 wind asymmetry during the minimum solar activity is not as larger as during the maximum solar 

316 activity. Also in Fig. [5] some stations located in the Northern hemisphere presented a larger 

317 correlation with VSS what can be attributable to the latitudinal dependence of the Sq current. 

318 Magnetic stations located at lower latitudes as VSS presented larger correlation (HON and SJG). 

319 Fig.|6]and|7]present the grafics of the diurnal and semidiurnal geomagnetic variations, respec- 

320 tively. These figures show the determination coefficient for the 12 chosen stations during the 

321 equinoxes and solstices between the years of 1999 and 2007. In the vertical axis, we present the 

322 determination coefficient and in the horizontal axis, the twelve chosen magnetic stations ordained 

323 in latitude from the southern to the northern hemisphere. This distribuion of stations helped to 

324 verify the latitudinal dependence of the diurnal variations. To validate their seasonal behaviour, 

325 we obtained the determination coefficient on March, June, September and December. In Fig. [6] 

326 and [7] March is displayed at top left, June at top right, September at bottom left and December 

327 at bottom right. We performed this analysis for the interval from high (1999) to low (2007) solar 

328 activity. This interval enabled the study of the solar activity dependence upon the Sq variations. 

329 Years of high solar activity are shown in black (corresponding to the mean value of the deter- 

330 mination coefficient for the years of 1999, 2000 and 2001), years of medium solar activity, in 

331 grey (corresponding to the mean value of the determination coefficient for the years of 2002, 

332 2003 and 2004) and years of low solar activity, in white (corresponding to the mean value of the 

333 determination coefficient for the years of 2005, 2006 and 2007). 

334 In Fig. [6] it was possible to identify a primary latitudinal dependence for the equinoxes and 

335 solstices. If we consider the simple dipole representation of the Earth's main magnetic field and 

336 the north-south symmetry of thermospheric wind and conductivity field between hemispheres in 

337 equinox conditions, we would find that the stations at the same latitude should have the same de- 

338 termination coefficient. However, the geomagnetic field is asymmetric so the mean con ductance 

339 for the northern and southern hemisphere will have different values even at equinox dSteningl 
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197 lh . This characteristic may explain the difference of response between some magnetic sta- 
tions located at similar magnetic latitudes but at different hemispheres. 

The disparity in the response of magnetic station between solstices is expected. These differ- 
ences between the d etermination coefficients could be explained by the following features of the 
Sq field reported by Matsushit a and Maedal ( 1965 ): 1) the total Sq current intensity is about 1.5 
times larger in the hemisphere's summer than in the winter and 2) the latitudinal position of the 
external current center is higher in summer than in winter. 

In this work, we also found that the 24-h diurnal variation presented differences in the determi- 
nation coefficient for magnetic stations located in approximately the same geomagnetic latitude 
but different longitudes. This mean that these stations were affected differently by the E-layer 
dynamo. At equinoxes, the determination coefficient was usually larger in years of higher solar 
activity, and at solstices it was practically larger in the year of lower solar activity. These re- 
sults may be explain ed by the variation of cond uctivity of the E-region which has a ionization 
density dependency (For bes and Lindzen , 19761) . Also, the total geomag netic field affe cts the 
conductivity through it dependence on electron and ions gyro-frequencies ( Steningi 1971 ). 

Another strong influence, very important to refer, is the intensity of the geomagnetic field. The 
total geomagnetic field is particularly high in the re gions of central Canada, Siberia and south 
of Australia and it is quite low near southern Brazil dCampbellL 1997 ). In Fig. [6] it was possible 
to observe that the magnetic stations, located at the northwest of Canada (CMO), north of the 
United States (BOU), China (BMT), Japan (KAK) and the central of Australia (ASP), where the 
geomagnetic field is particularly high, presented lower determination coefficient when compared 
to other stations. 

For the semi-diurnal variation, the latitudinal dependence was not as perceptible as in the diur- 
nal variation (see Fig.[6]and[7]i. Also, the equinoxes and solstices presented irregular distribution 
due to the solar activity. 

We expected lower correlation between VSS and CMO, because CMO is the station located at 
the highest latitude. The influence of the Sq currents is limited at high latitudes when compared 
to low and mid-latitudes. At high latitudes, the major influence comes from a large horizontal 
current that flows in the D and E regions of the auroral ionosphere, called the Auroral Electrojet. 
There are two characteristic of the auroral region that we must mention. First of all, the conduc- 
tivity in the auroral ionosphere is generally larger than at lower latitudes. Second, the horizontal 
electric field in the auroral ionosphere is also larger than at lower latitudes. Since the strength 
of the current flow is directly proportional to the vector product of the conductivity and the hor- 
izontal electric fiel d, the auroral e lectrojet currents are generally larger if compared to those at 
lower latitudes (see ISizoval (120021) for more details). 

Magnetic stations a nd observa t ories o n islands or coastlines are fully exposed to anomalous 
effects of the ocean (Schmucker, Il999allbh . The oceans, as well as the ionosphere, constitute 
electrical conductors and they are subject to tidal motions. As a result of the tid al flow of the 
ocean water, we might expect the effects of some kind of dynamo dParkinsorl 1983 ). The electric 
current system induced by the oceanic tides due to the drifting motion produces some kind of 
daily variation on magnetic stations located at oceanic islands and nearby shores. This shows 
that caution must be taken when interpreting the solar magnetic variations at a particular station, 
if the purpose is to obtain a world-wide analysis of the Sq field. 

The total variation measured on the ground c onsists of both external (ion ospheric current) and 
internal (induced Earth current) contributions (Forb es and Lindzenl fl976l) . The pattern of the 
conductivity of surfaces layers of the Earth will introduce a corresponding small scale pattern 
into the distribution of the induced currents. At certain stations, an accurate indication of the 
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average induc ed current sy stem is difficult to determine, consequently its effects can not be fully 
distinguished (Price, 1969). Howev er, the intensity of the external origin of the Sq field is about 
2.5 times that of the internal origin (Ma tsushita and Maedal 1 1965b . Therefore, in this work, we 
disregarded the influence of the internal contributions on the analysis of the diurnal geomagnetic 
variations. 



392 6. Conclusions 

393 A world-wide distribution of the harmonic components of solar diurnal variations (24 and 12-h 

394 periods) has been observed in latitude and longitude by using the gapped wavelet analysis and the 

395 wavelet cross-correlation technique. The objective of this paper was to study the characteristics 

396 of these variations at a Brazilian station as compared to the features from other magnetic stations 

397 to better understand the dynamics of the diurnal variations involved in the monitoring of the 

398 Earth's magnetic field. The main results in this analysis can be summarized as follows: 



399 1. The difference in the mean conductance due to the asymmetry of the geomagnetic field 

wo might be responsible for the difference in the latitudinal response between the magnetic 

401 stations, even at the equinoxes. 

402 2. We observed lower determination coefficients in the magnetic stations located on areas of 

403 the terrestrial globe where the geomagnetic field is particularly high. 

404 3. For the diurnal variation, the determination coefficient was usually larger in years of higher 

405 solar activity at the equinoxes and practically smaller in years of lower solar activity at 

406 solstices. 

407 4. The dependency of the E-region conductivity on the ionization density might explain the 

408 variation of the determination coefficient due to the solar activity. 

409 5. The variation of the E-region conductivity according to the intensity of the geomagnetic 

410 field might account for the differences of the determination coefficient among different sta- 

411 tions. 

412 6. The latitudinal dependence for the semi-diurnal variation was not as noticeable as for the 

413 diurnal variation, and also, it presented an irregular distribution due to the solar activity. 

414 7. In closing, this analysis showed that there were unusual and unequal contributions for each 

415 station according to their correlation to VSS. 

416 8. The analysis performed has shown a great spatial and temporal variability of the diurnal 

417 variation. These results suggest that the quiet days patterns and the physical processes 

418 involved still remain an important issue, that is far from being fully understood. 

419 The gapped wavelet analysis is an alternative technique to study data series that contain gaps 



420 of various sizes. Consequently, it can be used to exclude disturbed days and to study the char- 

421 acteristics of solar diurnal variations. Using a network of magnetic stations, we determined the 

422 contributions of the longitude and latitude and also of the total geomagnetic field. 
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a) Synthetic signal analysis 
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b) Synthetic signal analysis with addicional random gaps 
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Figure 2: Two spectrograms of synthetic signal consisting of a daily magnetic variation: (a) the synthetic signal and (b) 
the synthetic signal with addicional random gaps. From top to bottom, each spectrogram shows: the synthetic signal 
(signal with gaps), the continuous wavelet (gapped wavelet) spectrogram which satisfactorily resolves both in time and 
frequency and the phase change of the data. 
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d) Wavelet Phase 
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Figure 4: Each panel shows: (a) H component time series of VSS Station at June, 2007 with addicional gaps due to the 
disturbed days; (b) the local wavelet power spectrum using Morlet wavelet, logarithmic scaled, (c) the global wavelet 
spectrum and (d) the phase spectrum. 
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Figure 5: Modulus of the wavelet cross-correlation functions for each pair formed by VSS and one of the twelve chosen 
magnetic stations for June, 2007. 




■ set/HSA 

□ set/MSA 

□ set/LSA 




EYR AMS HER ASP HON SJG BMT BOU KAK CLF BEL CMO EYR AMS HER ASP H0N SJG BMT B0U KAK CLF BEL CM0 

i 100 - 



■ JUIVHSA _ go 

□ julYMSA 

□ jun/LSA 80 



■ dec/HSA 

□ dec/MSA 

□ dec/LSA 



70 















J3 






m 







EYR AMS HER ASP HON SJG BMT BOU KAK CLF BEL CMO EYR AMS HER ASP HON SJG BMT BOU KAK CLF BEL CMO 



Figure 6: Graphs of the diurnal geomagnetic variation during the equinoxes (March and September) and solstices (June 
and December) months, and the years from 1999 (high solar activity) to 2007 (low solar activity). The vertical axis shows 
the mean of determination coefficient and the horizontal shows the ABB code of the stations distributed by latitude. 
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Figure 7: Graphs of the semidiurnal geomagnetic variation during the equinoxes (March and September) and solstices 
(June and December) months, and the years from 1999 (high solar activity) to 2007 (low solar activity). The vertical 
axis shows the mean of determination coefficient and the horizontal shows the ABB code of the stations distributed by 
latitude. 
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